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ABSTRACT 

Most low-mass protostars form in clusters, in particular high-mass clusters; however, how low-mass 
stars form in high-mass clusters and what the mass distribution is, are still open questions both in 
our own Galaxy and elsewhere. To access the population of forming embedded low-mass protostars 
observationally, we propose to use molecular outflows as tracers. Because the outflow emission scales 
with mass, the effective contrast between low-mass protostars and their high-mass cousins is greatly 
lowered. In particular, maps of methanol emission at 338.4 GHz (J=7o-6o A+) in low-mass clusters 
illustrate that this transition is an excellent probe of the low-mass population. We here present a model 
of a forming cluster where methanol emission is assigned to every embedded low-mass protostar. The 
resulting model image of methanol emission is compared to recent ALMA observations toward a high- 
mass cluster and the similarity is striking: the toy model reproduces observations to better than a 
factor of two and suggests that approximately 50% of the total flux originates in low-mass outflows. 
Euture fine-tuning of the model will eventually make it a tool for interpreting the embedded low-mass 
population of distant regions within our own Galaxy and ultimately higher-redshift starburst galaxies, 
not just for methanol emission but also water and high-J CO. 

Subject headings: astrochemistry — ISM: jets and outflows — line: profiles — stars: formation — 
stars: winds, outflows 


1. INTRODUCTION 


Mo st stars form in high-mass clusters (|Lada & Lada 
2003), and the luminous star-forming regions seen m 
other galaxies are predominantly (embedded) high-mass 
clusters. In terms of luminosity, the cluster emission is 
dominated by the most massive O and B stars, and often 
global star-formation rates are infer red from their spec¬ 
tra a t optical and UV wavelengths (Kennicutt & Evans 
2012). Yet even the most massive clusters form very few 
high-mass stars compared to the total number of stars: 
the inferred star-formation rate thus depends on the mass 
distribution in the form of an initial mass function (IME) 
to acco unt for the contr ibution from the low-mass popu¬ 
lation (Hao et al.||2011). This population is notoriously 
difficult to access observationally even in G alactic high- 
mass clusters beyon d Orion ((i = 415pc,e.g., Hillenbrand 
& Hartmann|1998|) , in particular when this population is 
still forming and is deeply embedded such as in starburst 
galaxies. 

Nearby low-mass clusters, on the other hand, are read¬ 
ily studied and both the initial stellar mass and the core 
mass functions are well-constrained (e.g. Andre et al 


2010 Bastian et al.|[2MQ Offner et aL||2Q14[ ). However, 


it IS unclear if the formation conditions, and thereby the 
resulting core mass function, are the same in high-mass 
clusters where protostars are exposed to strong radia¬ 
tion fields from high-mass stars. This poses a problem 
when extrapolating conditions to extragalactic regions 


where individual embedded protostars are not resolved. 
Deeply embedded protostars, both low- and high-mass, 
are bright in molecular emission an d star-forming region s 
are now traced out to redshifts > 6 (Riechers et al.|2013), 
with molecular emission in local galaxies resolved into 


indiv i dual massive clust ers (e.g. Meier & Turner 2005 


2012 Meier et al. 2014). Thus, there is a need lor con- 


necting the nearby, well-resolved low-mass clusters with 
first the Galactic massive clusters, and then extrapolat¬ 
ing these to nearby extragalactic regions and eventually 
the higher-redshift starburst galaxies. 

One method for probing the embedded low-mass popu¬ 
lation is through outflow signatures, one of the key sign¬ 
posts of embedded star formation. The advantage of us¬ 
ing outflow emission is that (i) low-mass outflow emission 
scales with envelope mass (e.g. |Wu et al.|2QQ4[ |de VillieTsI 
et al.|2014| ), a relation which may also hold lor high-mass 
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protostars, thereby providing a lower contrast compared 
to the high-mass protostars in the region, and (ii) molec¬ 
ular outflows are uniquely associated with the embedded 
stages of star formation, both low- and high-mass. This 
combination makes outflows a potentially excellent probe 
of the low-mass end of the IME. 

While low-mass outflow properties (mass, outflow 
force) are typically measured based on spatially resolved 
low-J GO emission, CO also traces bulk cloud properties 
and is therefore not an ideal tracer of outflows and their 
parent population. Instead, molecules such as methanol 
are more suitable: the methanol abundance jumps by or¬ 
ders of magnitude in molecular outflows due to sputtering 
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(e.g. Bachiller et al. 1998) and in particular the mid-J 


transitions with K appear to be ide al from 


an empirical point of view (Kristensen et al. 2010). These 
transitions are not excited in the cloud where mostly low- 
J transitions emit K) and they are only 

moderately excited in hot cores where emission from 
higher-J transitions is more dominant (^up/te^lOO K)- 
The J=7k-^k transitions of methanol around 338 GHz 
are particularly well suited, both because they are good 
outflow tracers, and be cause instruments exist to map 
emission efficien tly (e.g. Kristensen et al.|2010 Torstens- 
son et al.||2011|). 

In this proof-of-concept study, we present a cluster 
model which simulates low-mass outflows and their emis¬ 
sion based on the known methanol emission from single 
low-mass protostars in the regions NGC1333 and Serpens 
Main. The outcome of the model is compared to observa¬ 
tions of an embedded high-mass cluster as a step toward 
calibrating the model before employing it for interpreta¬ 
tion of extragalactic observations, in particular by assess¬ 
ing what the high-mass outflow and hot-core contribution 
to the methanol emission is. The model is described in 
Sect. 2, and the results are presented in Sect. 3. The 
model is compared to observations in Sect. 4, where also 
the implications and future steps are discussed. 

2. MODEL 

The cluster model consists of two parts: (i) a template 
cluster with stellar spatial, age, and mass distributions 
and (ii) molecular emission assigned to each embedded 
protostar in the cluster. Below, these components are 
described separately. The codes for running the model 
are publicly availably An example of a cluster is shown 
in Fig. along with Hie flnal mass distribution. 

2.1. Cluster characteristics 

Spatial distribution. The spatial distribution of pro¬ 
tostars builds on the model by Adams et al. (2014) where 
the size of the cluster is i^max = /jS/q)^ and N 

is the number of stars. This relation works well with 
Ro=l pc, A^o=300 and 0 = 1/3 and these are the param¬ 
eters employed here. The radial probability distribution 
function is dP/dr=(3-p)/i?max (i*/Pmax)^~^ where the 
power-la w index p is s e t to 1 .5, again following the ex¬ 
ample of I Adams et al. (|2014). The angular distribution 
is assumed to be uniformly random. 

Higher-mass stars are often located closer to the bot¬ 
tom of the gravitati onal potential well, i.e. closer to the 


cluster center (e.g., Hillenbrand & Hartmann 


1998). To 


simulate this effect in the cluster model, the distance to 
the center is re-scaled by once the spatial dis¬ 

tribution is calculate d above, a method similar to the 


method employed by Moeckel & Bonnell (2009) where 
the index is 0.15. This small power-law scaling ensures 
that the most massive members of the cluster are located 
within 10-20% of the maximum cluster radius, whereas 
the lower-mass members are hardly affected. 

Age distribution. The age distribution follows that 


an estimated age of Myr (Evans et al. 

2009 

). The 

fraction of protostars in the Class I, H and 

HI stages as 


^ https://github.com/egstrom/cluster-in-a-box 
D01:lU.5281/zenodo.l3184 
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Figure 1. Spatial distribution of protostars in the cluster-in-a- 
box model. This particular cluster contains 3000 stars. The symbol 
sizes are scaled to the mass of each (proto)star. The inset shows the 
mass distribution after the low-mass population has been assigned 
ages and the Class O/I population assigned new envelope masses. 

well as the “ffat -spectrum” stage be tween Glass I and H 
~ where as the Glass 

(20141). Perseus is 


are taken from Evans et al. 


0 population is from badavoy et al. 


(20091 


one of the only low-mass clusters where such complete 
statistics are available. Where outffows only probe the 
embedded stages, the more evolved stages are revealed 
through, e.g., infrared obser vations where good statistics 
exist for the nearby clouds (Evans et al.||2009|). 

Mass distribution. The initial ma ss function (IME) 
follows that deflned in (|Chabrier|2003|) for young clusters 
and disk stars and is randomly sampled over the mass in¬ 
terval from 0.01-100 Mq (see inset in Eig. [B. Typically 
a cluster with 3000 stars contains ^10 hi^-mass stars 
(M > 10 Mq) and 1000 brown dwarfs (M < 0.05 Mq). 
Glass 0 cores typically have core masses t hat are three 
times higher than the flnal stellar mass (e.g. | Andre et al.| 
2010 and references therein), and Class I cores have to¬ 
tal masses twice as high as the flnal stellar mass. To 
account for this excess mass, the total mass of any core 
in these evolutionary stages is tripled or doubled, respec¬ 
tively. The masses of protostellar systems in the remain¬ 
ing evolutionary classes are typically close to their flnal 
mass but since molecular outflows from Class H and HI 
objects are rare and weak, these are not included in the 
model. 

In this model, the initial number of stars is specifled 
and the total mass of the cluster is subsequently calcu¬ 
lated. The stellar mass is from the IME (see above) and 
the protostellar masses are obtained by scaling the flnal 
stellar mass. The cluster mass is ^10^ Mq for a cluster 
with 3000 stars when all protostellar contributions are 
summed up. If the star forma tion efficiency is of the or¬ 
der of 10% ( iLada et al.||2010 ), i.e. the total stellar mass 
divided by the initial cloud niass, the initial cloud would 
have a mass of the order of ^10^ Mq which is typical of 
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observed Galactic molecular clouds. 

Outflow angular distribution. The outflow posi¬ 
tion angle is randomly set between 0° and 180°. Each 
outflow lobe is given a random distance from the proto¬ 
star to simulate both differences in inclination angles and 
dynamical ages, where the maximum separation corre¬ 
sponds to the maximum separation in the observed data 
(2x10^ AU; Table 0. 

2.2. Outflow emission characteristics 

The outflow template data are obtained from two low- 
mass star-forming regio ns, NGC1333 and Serpens Main 


(Kristensen et al 
Clerk Maxwell I'e 


2010|), both mapped with the James 
escope (JGMT) at 338 GHz. Three 


protostars are contained in the NGG1333 map (IRAS2A, 
4A, and 4B) and four in the map of Serpens (SMMl, 3, 
4, and S68N). The integrated intensity maps are shown 
in Fig. 

In this study, focus is placed on the strongest methanol 
line around 338 GHz, the 7o-6o A+ line at 338.409 GHz 
(Eup//cb= 65 K). The emission maps show empirically 
that this transition is an excellent outflow tracer, both 
because of the strong, high-vel ocity line wings and th e 
spatial distribution (Fig. j^and 


Kristensen 


er^[2^ 


This is in contrast to lower-excited lines, e.g. the J=2 
1 transitio ns around 96 GHz which trace the ambient 


Ltionj 

cloud (e.g. jOberg et al. 
transitions which typicaf 


2009), and the higher-excited 


tran sitions wnicn typicai r 
(e.g. |Menten et al. |l986 ). 
Observations of the J = 


y are associated with hot cores 


= 7o-6o A+ transition at 338.409 
GHz were obtained with the HARP instrument on the 
James Glerk Maxwell Telescope (JGMT) on Manna Kea, 
Hawaii. HARP is a 16-element array receiver which 
makes mapping observations very efficient. Both regions 
were mapped at an angular resolution of 15" over ar- 
eas of ^ 5^ x 5b The Serpens data are presented in 
Kristensen et al. (2010) where also the data reduction 
is described in detail. The NGG 1333 data were ob¬ 
tained from project M09AN05 (PI: Kristensen) and were 
reduced in a similar manner. The integrated intensities 
are converted from the T\ scale to flux by using a con¬ 
version factoij^of 15.625 Jy K 

Outflow emission from this transition correlates with 
envelope mass (Fig. [^. To test the strength of the cor¬ 
relation, a t test is performed and at the 90% confidence 
level there is a direct proportionality between envelope 
mass and intensity. Because of the low number of data 
points, the proportionality was not tested for higher lev¬ 
els of significance. The proportionality is S' = (78T33) 
Menv Jykms“^ where the uncertainty is Icr. 

To estimate the size of each outflow lobe in the model 
image, the sizes of the observed outflow lobes are first 
measured. This is done by measuring the area where the 
intensity is >10% of the peak, then deconvolving this 
area with the JGMT beam. Typical deconvolved radii 
are 1000-3000 AU, with a mean of 2000 AU. At a dis¬ 
tance of 3 kpc, this radius corresponds to 078 and so is 
resolvable with ALMA or the SMA, the latter in its very 
extended configuration. The emission distribution is as¬ 
sumed to be Gaussian for each lobe. The amplitude of 
the Gaussian is set such that the area matches the ob- 


Table 1 

Methanol observational parameters. 


Source 

-Menv^ 

(^o) 

ibol^ 

(.Lq) 

/ /UviB du dr2 
(K km s“^) 

(^pix) 

IRAS2A - East 

5.1 

35.7 

186 

23 

IRAS2A - West 

5.1 

35.7 

109 

27 

IRAS4A - North 

5.6 

9.1 

88 

16 

IRAS4A - South 

5.6 

9.1 

170 

21 

IRAS4B - Total 

3.0 

4.4 

34 

13 

SMM4 - South 

6.8 

6.2 

76 

21 

SMM4 - West 

6.8 

6.2 

59 

15 

SMM3 - Total 

10.4 

16.6 

22 

12 

S68N - South 

II.7 

II.7 

148 

21 

S68N - North 

II.7 

II.7 

136 

30 

^ From iKristensen 

et al.|(|2012|, except for S68N which is from |Kris-| 

|tensen et al.|(|2010||. The parameters for the Serpens sources 

have been 

scaled to a distance of 415 pc 

(|Dzib et al.|2010||. 


^ Number of 7 '.' 6 X 

7/5 pixels. 




served flux. All observed parameters are listed in Table 

m 


2.3. Outflow emission assignment 

Each low-mass protostellar object is assigned outflow 
emission based on a scaling of the observed emission with 
modeled envelope mass; that is, for a given model enve¬ 
lope mass, each simulated outflow lobe is given a total 
integrated flux corresponding to the scaling of the ob¬ 
served fluxes. 

The outflow activity of Glass I sources is typically 
about a n order of magnitude lower than for Glass 0 
sources (Bontemps et al. 1996). This lower activity is 
included in the model by scaling down the emission by a 
factor of ten for each Glass I outflow. Finally, the map 
is convolved with the desired beam. Note that no radia¬ 
tive transfer is done explicitly, and that all emission is 
impli citly assumed to be optically thin (Bachiller et al. 


3. RESULTS 

For a cluster with 3000 protostars and an age distri¬ 
bution of the low-mass population corresponding to that 
of Perseus, i.e. 1 Myr, the resulting methanol emission 
map is shown in Fig. The image has been convolved 
with a beam of 075 and the cluster is assumed to be at 
a distance of 3 kpc. 

Almost all emission is located in a central condensation 
where individual outflow lobes are marginally resolved 
as little clumps (the FWHM of each outflow lobe is 078; 
see above). The maximum separation between outflow 
lobes is 2x10^ AU, or ^ 7" at a distance of 3 kpc; since 
each outflow lobe is treated separately as 2D Gaussian 
functions, two lobes will only be identified as originating 
in the same protostar if the lobes are close, ^2". For a 
larger beam, such as ^ 2" which is easily achievable with 
the SubMillimeter Array (SMA), the outflows would have 
been unresolved. 

The peak intensity is 0.4 Jy km s“^ beam“^ in the 
model image. Figure shows the radial flux profile. The 
profile peaks 3" from the center and remains above 0.5 
Jykms“^ to radii >10" (Fig. dashed line). 


^ http://docs.jach.hawaii.edu/JCMT/HET/GUIDE/het_guide/ 


4. DISCUSSION 
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Figure 2. Methanol emission toward NGC 1333 (top) and Serpens Main (bottom) obtained at the JCMT. The gray scale is for flux 
and the contours approximately outline the regions over which outflow emission was integrated per outflow lobe, with red contours for 
redshifted and blue for blueshifted lobes. In the cases where two lobes are denoted (North/South, East/West), the protostar is located in 
between them (the protostars are marked with stars). The lobes labeled “total” include spatially unresolved emission from both the red- 
and blue-shifted lobes. The gray lines show the extent of the maps. 


4.1. Comparison to ALMA observations of 

IRASI6547-4247 

The high-mass star-forming regi on IRAS 16547-424 7 
was recently observed with ALMA (Higuchi et al.|2015). 
The region is located at a distance or 2.9 kpc and was ob- 
served at an angular resolution of ^ 0'.'5 (1500 AU). The 
data were downloaded from the Science Archive. The 
pipeline-calibrated visibilities were were first continuum- 
subtracted using line-free channels in the uv spectrum 
and then re-imaged; two changes were made compared 
to the pipelined image presented in Higuchi et al. (2015): 
first, the visibilities were averaged over a spectral re- 
gion corresponding to velocities from -44 to -20 kms“^ 
prior to imaging; second, the clean mask was iteratively 
placed so that cleaning started at the peak intensity and 
expanded from there. Both were done to increase im¬ 
age fidelity in the final restored image, with the added 


bonus of a decrease in noise level of a factor of 2 to 0.07 
Jy kms“^ beam“^. The ALMA data are shown in Fig. 
^ where they are not corrected for the shape of the pri¬ 
mary beam. For the more quantitative analysis, see be¬ 
low, primary-beam-corrected data are used. The same 
primary beam is taken into account when analyzing the 
simulated data. 

The emission is spatially resolved into a central conden¬ 
sation, and a larger number of sub-structures. The cen¬ 
tral source drives a large-scale CO outflow and part of the 
methan ol emission is clearly associated to the high-mass 
outflow ( Higuchi et al.|2QT5 ). The more wide-spread dif¬ 


fuse methanol emission is likely connected to outflow 
activity based on the large line-width (^10 kms“^). 


(Brown & Bolina 

|2007|), and such a 

^ a radius of <1" (1 

disschop et al.||2007) 


Thus 
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Figure 3. Methanol emission versus envelope mass, Mgnv, for the 
sources in NGC1333 and Serpens. The error bars are for the 30% 
calibration uncertainty, and the line shows the best-fit proportion¬ 
ality. 


it appears likely that the spatially extended emission is 
caused by other , perhaps low-mass, outflow sources as 
also argued by (Higuchi et ah 2015). Below, the stan¬ 
dard model results are compared to these observations 
without fine-tuning of the input parameters. 

To compare the model and observations in a quanti¬ 
tative manner, both images are radially averaged from 
the center of the cluster including only fluxes >3cr. The 
same observationally set rms cut-off was also applied to 
the model image. The resulting radial profiles are shown 
in Fig. [dj The signal at distances >8" may be diluted 
by false-positives at the 0.1-0.2 Jykms“^ level. At dis¬ 
tances below ^2", fluxes are dominated by the high-mass 
source (outflow and hot core), and at larger distances the 
model profile is 30-50% weaker than the observed profile. 
If the rms cut-off is not applied to the model image, the 
radial profile overlaps with observations to >10" without 
a scaling. 

Assuming that the emission at distances > 2" is domi¬ 
nated by the low-mass population, the total flux fraction 
of the low-mass outflows from the entire cluster can be 
estimated: approximately 50% of the cluster flux is com¬ 
ing from the low-mass outflows, the rest is from a com¬ 
bination of the molecular hot core and the high-mass 
outflow(s). Disentangling these contributions is beyond 
the scope of this work. 

To obtain the factor of two increase in the model fluxes, 
several possibilities exist: 


• Age. If the cluster is younger than 1 Myr, the rel¬ 
ative number of Class 0 and I sources increase and 
so will the number low-mass outflows. The model 
presented in Fig. [^has ^150 and 400 Class 0 and 
I sources, respectively, spread over a radius of 150" 
(but of course centrally peaked). Doubling these 
numbers would lead to a flux doubling as well. 

• IMF. If the IMF is shifted toward higher masses, 
which could be the case if the higher-mass stars 
form before the lower-mass ones, then more 
methanol emission would be generated per object, 
because each object would be relatively more mas¬ 
sive than in the standard model. 


• Size. The observed cluster may be forming more 
than 3000 stars in which case more Class 0 and I 






Offset (arcsec) 



Figure 4. Top: Integrated intensity map of methanol emission at 
338.4 GHz. The image was convolved with a O'.'S beam. The color 
scale is shown on the right. Center: AL MA image of the 338. 4 GHz 
emission toward IRAS 16547-4247 from|Higuchi et al.|(|2015|). Bot¬ 
tom: Radially averaged flux from the center ot tlie cluster (blue), 
multiplied by 2. Also shown is the radially averaged flux profile of 
the ALMA image (green). Only fluxes >3(7 are included in these 
profiles. The dashed line shows the radial model profile without 
(7-flltering and without scaling. 

protostars would need to be included in the model. 


• Spatial distribution. The simulated Class 0 and I 
protostars are spread randomly over the entire clus¬ 
ter: if the stars form closer to the center, the pro¬ 
tostars would be more concentrated and the flux 
would increase. Indeed, the median distance for a 

























6 


Kristensen & Bergin 


Class 0 protostar to the cluster center is 0.4pc, or 
25". 

All of these possibilities will be explored quantitatively 
in the future. However, given that the model is not tuned 
to reproduce the observations such an overlap is excellent 
and lends credibility to the model predictions even when 
using a standard IMF. 

4.2. Future prospects and outlook 

This proof-of-concept study demonstrates one possi¬ 
ble method for gaining access to the embedded low-mass 
protostellar content of a forming cluster, namely through 
molecular outflow emission from these low-mass proto¬ 
stars, and thereby potentially constraining the initial 
(core) mass function. Because outflow emission is di¬ 
rectly proportional to the envelope mass, the contrast 
between low- and high-mass protostars is effectively low¬ 
ered by several orders of magnitude compared to the 
bolometric luminosity. 

The success of the model in reproducing part of the 
morphology of a single observed high-mass embedded 
cluster paves the way for further fine-tuning of the model, 
and for including the contribution from high-mass out¬ 
flows and hot cores explicitly in the model. Two avenues 
will be pursued in the near future: (i) parameter space 
will be explored, in particular with respect to constrain¬ 
ing the initial (core) mass function of low-mass embedded 
protostars, and (ii) other species will be included. As 
for the latter, Herschel observations demonstrate that 
H 2 O and high-J CO (J>16) trace shocks and outflows 
uniquely, and these are the next species and transitions 
to be included. Although no new observations of H 2 O 
will be done for the foreseeable future, Herschel already 
observed hundreds of Galactic star-forming regions, both 
low and high mass, thus providing enough data to cali¬ 
brate the model. Furthermore, ALMA will easily be able 
to observe H 2 O emission in redshifted galaxies, and mod¬ 
els of this type carry the promise of providing constraints 
on how low-mass stars have formed and how the IMF has 
changed throughout cosmological times. 

We thank the referee for thoughtful and constructive 
comments that helped to improve this manuscript. 
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